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Abstract. 

We define infinitesimal contact transformations for ordinary difference 
schemes as transformations that depend on K + 1 lattice points (K > 1) and 
can be integrated to form a local or global Lie group. We then prove that such 
contact transformations do not exist. 
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1. Introduction 

Let us consider a system of ordinary or partial differential equations involving the 
independent variables x € TZ P and the dependent variables y £ lZ q . Three types of 
transformation of variables play important roles in the theory of differential equations. 

The first are point transformations, in which the new variables (x, y) depend only 
on the old variables (x,y). The methods of Lie point symmetry groups have been 
adapted to make them applicable to difference equations ( and differential - difference 
equations [TJHISHE]- 

The second are generalized transformations for which (i, y) depend not only on 
x and y but also on derivatives of y with respect to x, like y a , Xi , ya,xix k , ••• ■ These 
are important in the theory of integrable partial differential equations. They too have 
been adapted to difference equations and other equations involving discrete variables 
(or variables on lattices) [51 M IT3]. 

Intermediate between point and generalized transformations are contact ones. As 
symmetry transformations they provide non trivial results only for scalar equations [5] 
For them we have 



(no second derivatives). 

The Lie algebra of contact transformations is closed under commutation, i.e. 
it does not introduce higher derivatives. This makes it possible to integrate the 
Lie algebra to a Lie group of transformations of the form (|T|) and ©. The 
inverse transformations have the same form. For some differential equations contact 
transformations provide analytical solutions, or at least a reduction of the equation, 
not obtainable by using point transformations [4"I [Tl"irB] . 

Contact transformations have, so far, not been adapted to difference equations. 
In this article we will just consider the case of ordinary difference equations. The 
purpose of this article is to analyze the situation in this case and show that contact 
transformations for ordinary difference equations do not exist. 

In Section 2 we will review the theory of contact transformations in the case of 
ordinary differential equations while in Section 3 we review the theory of Lie point 
symmetries for difference schemes. Section 4 is then dedicated to proving that contact 
transformations for difference equations do not exist. 

2. Point and contact transformations for ordinary differential equations 

Let us consider an ODE of order n 



V = tt(x,y,y x ), 



x = A(x,y,y x ) 



(1) 



and a contact condition assures that 



Vx = $(x,y,Vx) 



(2) 



yW =F(x,y,y',...,y^) 



(3) 



and the n-th prologation of an element of its symmetry algebra 



The coefficients in the prolongation are calculated recursively 




(5) 
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where D x is the total derivative operator. The symmetry condition is 

w in) X{y (n)_ F )\ y{n)=F = Q. ( 6 ) 

For point transformations we have (by definition) 

f = f(x,t) <f> = <Kx,y) (7) 

and hence eq. ([5]) implies c/)( n ' = ^"'(s, y, y', . . . , y( n '). A Lie algebra of point 
transformations can be integrated to a Lie group. Indeed we can integrate the vector 
field X 

§ = «2,y), % = (8) 

The transformation formulas for the derivatives y^j can be obtained (for any /c) 
either by the chain rule, once y(x) is known, or by integrating pr^ k 'X. 

Now let us consider generalized symmetries containing higher order derivatives 
upto order n. In this case, instead of ([7]), we postulate 

Z = Z(x,y,y',...,y<ri), <t> = <t>(x,y,y',...,y^) (9) 

To be able to integrate the Lie algebra of such transformations to a Lie group, we 
need the coefficients in the n-th prolongation (for k = 1, 2, n) to depend only 
on x, y, y', y^ but not on y( n +^) , j > 1. H. Stephani [10] provides an elegant proof 
of the fact that this condition cannot be satisfied for n > 1. 
For n = 1 contact transformations do exist if 

</y = (10) 

and an ODE can be invariant under a group of contact transformations. The elements 
of Lie algebra have the form 

X = ^(x, y, y')d x + <f>{x, y, y')d v + <f>W (x, y, y')d y , (11) 

The contact condition (JTDJ) eliminates y" from (fP^ and assures that c/>W for k > 1 
is given by (O . Thus first order contact transformations of the form (fTTj) form a Lie 
algebra that can be integrated to give a Lie group of contact transformations. 



3. Ordinary difference schemes and their point transformation 

A difference equation on a fixed non-transforming lattice has very few continuous 
symmetries. One way of making a Lie symmetry approach to difference equations 
fruitful is to consider difference equations on flexible lattices that themselves transform 
under the transformations [3[2J[5]. Thus instead of difference equations we will deal 
with difference schemes. 

Let us consider the case of one independent variable x and one dependent one y. 
The independent variable x is sampled at several points Xk and y(x) is evaluated at 
the same points yt = yk(xk)- An ordinary difference scheme (OAS) consists of two 
relations between the K point Xk,yk- 

E a ({x k }, {y k }, k = n + M, n + M + 1, . . . , n + N) = (12) 

a =1,2 K = N-M+1, n,M,NeZ, N > M 

In the continuous limit ajj+i — Xj — foj+i — > one of the two equations (|12p should go 
into an ODE of order K — 1 (or less), the other into an identity (like = 0). If the 
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values Xk , Uk ar e given in K — 1 neighbouring points then we must be able to calculate 
their values in the next point. The corresponding independence condition is 

J {El ' E2) ^ 0, Vn (13) 

0{X n+N , y n +N) 

Thus the OAS (JT2J) determines both the lattice and the solution of the difference 
equation. The general solution of (IT21) can be written as 

x„=x(n,ci,c 2 ,...,c 2 (K-i)) (14) 
Vn = y(n,ci,c 2 , . . . ,c 2 (k-i)) = 4> n {x n ,c 1 ,c 2 , ■ ■ . ,c 2 (k-1)) 

Point symmetries of the OAS (| 12[) can be determined algorithmically. We require 
that the elements of the Lie point symmetry algebra of ([12]) have the form 

X = £, n dxn + </>n9 y „ ( 15 ) 

with £„ = £(x n ,y n ),(f) n = <j)(x n ,y„) and the prolongation of the vector field (p~5)l 
should have the form 

n+jV 

prW X = i^d Xk + ^dyj (16) 

k=n+M 

where the sum is over all points figuring in the OAS (|12j) . The determining equations 
are obtained from the conditions 

pr^XE a \ El=0 , E2=0 = a =1,2. (17) 

Since each function and fa (for point transformations) only depends on one point 
{xk,Vk) equation (|17|) actually splits into a set of several functional equations for the 
functions and fa. 

If we fix the value of n then the K points in (fl"2"T) form a "stencil" [2] • Within a 
stencil we can choose different coordinates, more appropriate for taking a continuous 
limit. Let us consider a stencil with the points {xn+k—i, Vn+k-i, 1 _ k < if} for some 
fixed n. An alternative set of coordinates on the stencil is 1121: 



with 



r f 1 ) ( 2 ) ( 3 ) ( 4 ) {K-l) . . . /i Q s 



fl-n+k — %n+k ~ ^n+fc-l 

P% = (19) 

(2) _ ^n+2 Pn+1 



Pn+2 



= 2- 

^n+2 

(if-2) _ (if-2) 
(K-l) _ (K .s Pn+K-1 Pn+K-2 
Pn+K-1 - l A ^ 

^n+K— 1 

In the continuous limit we have 



If we transform the "discrete" prolongation (116)) of the vector field (1151) to the new 
variables (fl8]l . we obtain 

AT AT 

prX = £ n S Xn + <M„„ + £ « (fc) &„ + * + £ ( 21 ) 

fe=i fc=i 
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The general formulas for the coefficients and (j)^ are 

where the A T is the total difference operator 

A T F(a5 n ,2/n,p£ji,p£ 2 j. 2 , • • •) = 7— { ^(Zn+l, yn+l,pi+ 2 ,pi+3 » • • •) 

- ^Sn.i/n.pSi.pSa. •••)}■ ( 23 ) 
4. Do contact transformations for OAS exist? 

Let us again consider the OAS (IT2"j) . The question that we pose is the following. Is it 
possible to define contact transformations that leave the solution set of the OAS (fT2"| 
invariant? We shall impose the following restrictions to define contact transformations 
in the discrete case. 

Definition 4.1 A contact transformation for the OAS U2\) satisfies: 

(i) The vector fields forming the Lie algebra L c of contact and point transformations 
should have the form M5\) with at least one of the coefficients and (f>„ in one 
of the fields X depending on K + 1 points. 

Cn = £,n(x n ,y n ,x n+ i,y n+1: . . . ,x n+K: y n+K ), K EZ + 

4*n 4*n (*£n 3 Un i •En+l 3 Vn+1 3 • ■ • 3 •En+K 3 yn+K ) (24) 

(U) The Lie algebra should be integrable to a Lie group. This implies that all 
coefficients h22]) in the K — th prolongation of X h21)) should depend only on 
the points (x n , y n , . . . , x n +K, Uu+k) and not on any further ones. 

(Hi) In the continuous limit the algebra and the group of contact transformations of 
the OAS il8\) should reduce to the corresponding Lie algebra and Lie group of 
contact transformations of the corresponding ODE. 

Definition 4.2 A contact transformation for a difference scheme is said to be of order 
£ if it involves £ + 1 points of the lattice. 

Let us now prove the following theorem: 

Theorem 4.1 Contact transformations of order 1 of OAS do not exist. 
Proof 

Let us consider the first prolongation of a vector field depending on two points: 
pr^X = i n 3 Xn + 4> n d yn + n^d hn+1 + (25) 



We have 



Cn in^^ni Vni 2/ra+l) — £ra(*£nj Vni ~t~ ^n+ls Vn ~t~ ^n+lPn+l) 

<t>n = 4>n{x n ,yn,X n +l, Dn+l) = 4>n(%n,yn,X n + h n +i,y n + hn+lPf$§) 

,(1) _ 4>n+l - <K (1) gn+1 ~ jn (e} ^ 

Vn+1 — 7 Pn+1 T 3 \ z >) 

lin+1 n n +i 
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£n+l = Vn+li X n + 2 , Vn+Z) 

= £n+l{Xn + h n +l,y n + h n +lPn+li%n + h n +l + h n+2 ,y n + 

i ft. i i u (1) , h n+2 (2) n 

+ (hn+i + /l„ +2 )K+l + — g- Pn+2l) 

0ri+l = 0n+l (x n +i,y n +i,X n + 2 , Vn+2) (28) 

= 0n+l(^n + h n +i,y n + hn+lPn+D x n + h n+ i + h n +2,Vn + 

i ft. i 7 u (1) . (2) n 

+ (h n+1 + h n +2)\p n + 1 + — Y~ PnU\) 

As one can see in ([28]) both k^ 1 ' and depend on (/i„ +2 ,pi+2)- The question is 

whether this dependence can cancel out by proper choice of the symmetry generator 
form. We have: 



d&n+i _ h n+ 2(h n +i + h n+ 2) 



' J n+1 d£n+l 



dyn+2 Pn+1 dy,, 



Similarly dn^ /dh n+2 = and d^li/ dy n+2 = imply 



din d<j) n 



0, (29) 

. / ' i .-!_ I (III,,—-) ('II,.— -1 

which provide a contact condition similar to the continuous one (1101) . However in this 
case we have to require the further condition 

<9k (1) _ h n+2 (h n+1 + h n+2 ) <9£„ +1 _ 

9p n %~ ^ 8y n+2 -° W 

and hence d£ n +i / dy n + 2 = and consequently from (|29p d4> n +i/ dy n + 2 — 0. In turn 
these conditions imply 

-^ = ^ = 0. (31) 
oy n+ i ay n +i 

) 

= 0, (32) 
ax n +i ax n+ i 

but (|31[) and (|52"j) mean that ([23)1 corresponds to a point transformation, rather than 
a contact one. □ 
By a similar calculation one can prove that if £ and <fi depend on K points then 
the K — th prolongation will depend on one more point and hence one cannot obtain 
an integrable (closed) system of contact transformations for an OAS. To sum up we 
state: 

Theorem 4.2 Contact transformation of any order K > 1 for OAS do not exist. 

If we consider discrete equations on a fixed (nontransforming) lattice then most, 
or all, point symmetries are absent. In this case the variable x n — hn + xg with h and 
xq fixed numbers and thus the infinitesimal generator is just of evolutionary form, i.e. 

X = <j> n (y n , y n +i)dy n + (f> n +i(y n +i,y n +2)dy n+1 = <p n d yri + ^ 1 | 1 9^(33) 

Then (|33|) is a special case of (f2Tj) with and K,( k > absent. Thus the nontransforming 
lattice is a special case of the previous theorem. 
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5. Conclusions 

Theorem 14.21 amounts to a "no-go theorem" . It states that for ordinary difference 
schemes contact transformations, as defined in Definition 14.11 do not exist. This 
leaves open the possibility that generalized symmetries might exist that in the 
continuous limit go into contact symmetries of an ODE. For comparison we recall 
that when considering integrable differential-difference equations on fixed (non- 
transforming) lattices the following situation occurs. An infinite Lie algebra of 
generalized symmetries exists of which a small subset "contracts" to point symmetries 
in the continuous limit |3J. Lie point symmetries of the corresponding differential 
equation are recovered in this manner. 
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